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Abstract. 

Quantum lithography promises, in principle, unlimited feature resolution, 
independent of wavelength. However, in the literature at least two different theoretical 
descriptions of quantum lithography exist. They differ in to which extent they predict 
that the photons retain spatial correlation from generation to the absorption, and 
while both predict the same feature size, they differ vastly in predicting how efficiently 
a quantum lithographic pattern can be exposed. 

Until recently, essentially all experiments reported have been performed in such 
a way that it is difficult to distinguish between the two theoretical explanations. 
However, last year an experiment was performed which gives different outcomes for 
the two theories. We comment on the experiment and show that the model that fits 
the data unfortunately indicates that the trade-off between resolution and efficiency 
in quantum lithography is very unfavourable. 
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1. Introduction 

For a long time it was believed that the resolution of optical systems was limited 
by the so called Rayleigh criterion [1], which roughly states that the resolution of 
the system is diffraction limited by the half of the wavelength A of the used light. 
Hence, it was believed that only by going to shorter and shorter wavelengths one could 
improve the resolution substantially. However, going from the visibile to UV-light, 
x-rays, or gamma radiation means using light with such a high photon energy that 
the radiation may damage or destroy the object one wants to expose or investigate. 
Furthermore, it is difficult to build optical systems, since ordinary lenses are no longer 
transparent in these wavelengths regions, and even gases start to absorb the light quite 
substantially. Moreover, building mirrors are also challenging due to absorption and very 
small refractive index differences between different materials. This makes it complicated 
and expensive to build Rayleigh limited optical systems for wavelengths shorter than 
about 200 nm. 

There have been proposals for "classical" systems which achieve higher resolution 
than the Rayleigh limit [21 [3l H], but such lithographic schemes suffer from visibility 
problems and they cannot simultaneously generate a sub-Rayleigh feature size and a 
sub-Rayleigh feature separation, at least not until new kinds of detectors are developed 

It therefore came as a relief when Boto et al showed that entanglement could be 
used to beat the diffraction limit given by the Rayleigh criterion [8]. They showed 
that the use of entangled photons would allow one to surpass the diffraction limit 
by a factor of compared to the Rayleigh criterion. Lithographic systems would 
potentially profit substantially from the idea of Boto et al by allowing the generation 
of features arbitrarily smaller than the Rayleigh limit [9l[T0l[IIl[l2l[T3l[Tl[l5l[l6]. To 
make these systems work in practice, there is still a plethora of problems to solve. For 
example, one has to find a material which changes its behaviour only when absorbing N 
photons [21 [171 [iHj , and one needs also to optimise the absorption rate of this material 
[19j. Unfortunately, there are limits to the extent this is feasible [20]. Additional 
difficulties involve the rather peculiar, multi-photon, entangled photon states that need 
to be generated [211 122]. Alternatively, other ways to generate indistinguishable A^- 
tuples of photons have to be developed [23]. There have been investigations for the 
case of multiple frequencies or in the broadband limit [211 EH]. For some schemes even 
photon losses may become an issue [26]. Furthermore one has to think about how to 
generate arbitrary patterns [TOl [HI [121 [27]. An open question that we address in this 
paper is if systems built on this principle also would work efficiently in practice. The 
main aspect is whether or not the time needed to irradiate a lithographic film would be 
short enough to make this method efficient. 

We start our paper in section [2] by describing two different models for most of the 
experiments hitherto reported in the field of quantum lithography, but where one of 
the explanations predicts a favourable time scaling whereas the other one does not. We 
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Figure 1. Setup for quantum imaging or lithography for N — 2. Each of the in total 
S" + 1 coordinate points on the screen to the right consists of a detector with width 
h which detects two-photon absorption. The double slit on the left has a distance d 
between the slits and is at a distance R to the screen. 



then discuss an experiment in section |3] which would give different outcomes for the two 
explanations and which thereby could decide which of the explanations is correct (or at 
least applicable to this kind of experiment) and how the time for imaging or radiation 
scales with the number of photons N per state and the number of image "pixels" S. 
The outcome of an experiment recently performed by Peeters et al [28] allows us then 
to deduce which of the explanations seems to be applicable, and the result suggests that 
the time-scaling behaviour of such a quantum lithography scheme is, unfortunately, very 
unfavourable. In section H] we relate our paper to other papers in the field before we 
conclude our results in section [5l We have also two appendices, one where we give a 
thorough mathematical analysis and derivation of our model and one where we look in 
more detail into the original proposal of Boto et al . 

2. Different theoretical explanations 

Imagine a setup as depicted in figure [H A light beam of wavelength A is incident on a 
double slit with slit distance d':^ \. Due to diffraction at the slit, the light will produce 
an interference pattern on a screen at a distance d behind the slit. Assume that our 
screen now consists of S + l detectors, each of the same width b. Below, we have assumed 
that S is even (representing an odd number of detectors). However, the conclusions we 
draw are equally valid if an even number of detectors are used. The detectors cover 
the whole length of the screen and b is small enough to clearly be able to resolve the 
interference pattern, i.e., the distance between two intensity peaks in the detector plane 
is many times larger than b. A photon impinging on the double slit can either pass 
through the upper or through the lower slit, i.e., the state at the slit will be 




(1) 
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Since we have decided to remain ignorant as to which path the photon took through the 
double sht, but only detect the photon with the detectors at the screen, we will there 
have the state 



t=—— 
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where M\ = a/2(5' + 1), the summation is over all detector modes s, k is the wave vector 
of the incident light, r (s) = y/B^T~sW ^ R (^1 + and Ar (s) = dsin6' ^ d6' = 
dbs/R. (For the normalisation in expression ([2]) we have for simplicity assumed that the 
pattern falling on the 5 + 1 detectors contains exactly an integer number of fringes.) 
The ensuing measurement probability at the detection plane is 

, , 2 / ^^^^ \ / o \ 

Pl,,c(s) OC COS (^^)- (3) 

The distance between two peaks in the measured interference pattern is given by XoR/d, 
where Aq is the used wavelength. In the equation above we have neglected the envelope 
function induced by the exact shape of each slit. Instead we have assumed that the 
photons have equal probability amplitude to hit any detector. This is a simplification 
justified by the fact that it does not alter our main results or conclusion. 

If we look at the quantum lithography scheme the setup is slightly changed. We 
are not any longer interested in single photons but in "bunches" of photons. For 
most of our paper we restrict ourselves for practical reasons to = 2 photons, but 
our whole argumentations will be valid for any number A^. To show that one can get 
a higher resolution (i.e., a shorter distance between peaks) in this setup, one wants to 
assure that two indistinguishable photons pass through the same slit, i.e., that one has 
a two-photon, NOON state at the slits [8] 

l^2,si) = ^ (|2)„ ® |0),„^,, + |0)^pp^, ® . (4) 

This can be implemented, e.g., by placing a short type-I crystal in front of the double 
slit and pumping this crystal so that it produces pairs of degenerate photons which leave 
the crystal collinearly. One also has to assure that the pump beam is broader than the 
distance between the slits and that the two slits are centered on the setup's optical axis. 
Note that the state at the screen is isomorphic to a Bell state, that is it has maximal 
transverse entanglement of its photon number. Also note that, since the slits can be 
viewed as the source illuminating the screen, an analysis of how this state is produced 



is not needed. Indeed, in Appendix A it turns out that even long crystals, narrow 



beams and a long distance between crystal and shts can produce |^E'2,si)- It suffices 
to conclude that the state |^'2,si) is the optimal 2-photon state for a demonstration 
of 2-photon interference fringes, giving unity two-photon interference visibility and 
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interference fringe doubling as compared to the classical Young's experiment at the 
same wavelength of light. 

The interesting question is what happens with the state \ '^2,si) after that the photons 
have passed the double slit. One could use the concept of photonic de Broglie- waves 
|29] to explain the effect of quantum imaging or lithography. Boto et al [8] analysed the 
probability for the two-photon state to arrive at a specific point on the screen depending 
on the intensity / of the incoming light as follows: "For two-photon absorption with 
entangled photon pairs, the absorption cross section scales as I... If the optical system 
is aligned properly, the probability of the first photon arriving in a small absorptive 
volume of space time is proportional to /. However, the remaining — 1 photons are 
constrained to arrive at the same place at the same time, and so each of their arrival 
probabilities is a constant, independent of /." 

Transferring this idea to the experiment in figure [H would give the state 



*~ 2 



where Mb = ^4:{S+1). At the screen, the distance between two peaks XqR/ {2d) is only 
half as long as in ([3]) and one could treat the two photons as a photonic de Broghe-wave. 
The implication of this state is that the two photons can hit essentially any detector 
(except those at the interference anti-nodes) but that both photons will hit the same 
detector. 

Steuernagel [30], however, opposed that description by stating that "... it is not 
true that the first arriving photon greatly constrains the arrival location of the following 
ones ... Very few photons will be absorbed in one point since they typically arrive far 
apart." In a mathematical description his argumentation can be given the following 
description: 



2 2 



where 



^st.c) = ^ E E {e^'^^-^^)^-(*)latat + e''=[-(^)+-Wlatat} |0) , (6) 



-MW^flll + ^U^. (7) 



is the distance from the upper and lower slit, respectively, and A/st = V2{S + 1). 
Although one might argue that our derivation of ([5]) and ([6]) is somewhat short-handed 
or hand waving, a complete quantum mechanical analysis of the state from the pump 



beam to the screen leads to the same result. Such an analysis is given in Appendix A 



Also Tsang opposed to the assertions of Boto et al . He based his conclusion 
on a rather involved analysis, and found that the multi-photon absorption-rate should 
actually be lower for spatially entangled states, such as NOON-states (of which the state 
|\l/2,si) is a A^ = 2 realisation), than for classical states [31]. However, in the following 
we will follow Steuernagel's somewhat simpler model. 
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The difference between the descriptions is that the photons retain a spatial 
correlation ("stick together") in the description of Boto et al , and hence always 
will arrive at the same position. In Steuernagel's model the photons will propagate 
independently of each other after the slit, and can hence arrive at different detectors. 
When looking at the case where both photons arrive at the same detector (that is, s = t, 
see below) the models coincide when it comes to the spatial behaviour and thus, by only 
looking at two-photon, same position events, one will se identical interference patterns, 
both compressed by a factor of 1/2 as compared to the single photon interference, 
described in ([2]) and ([3]). 

Although predicting the same two-photon, same-postition absorption pattern, the 
two descriptions in ([5]) and (|6]) will give a quite different scaling behaviour for how 
efficient a two-photon exposure would be. In (j5]) both photons in every two-photon state 
passing through the slits arrive at the same detector, and going from 2 to photons 
in (jl]), the exposure efficiency will be invariant. (However, the absorption efficiency of 
any detector or material will in general be considerably lower, but this effect will be the 
same in both cases). Changing the number of detectors S will only linearly increase the 
time one has to expose to get a certain number of A^-photon events per detector. Hence, 
the "exposure time" of a quantum lithographic image will scale as t-q{S,N) oc S. 

Equation ([6]), however, assumes that the arrival of the first photon at a certain 
detector does not in any way determine where the other photon will arrive and therefore, 
in the two-photon case, the exposure time rst(5', 2) oc 5*^. The photons in an A^ photon 
state will also hit the detectors in a random fashion, but only if all photons arrive 
at the same detector will the detector fire. Hence, the exposure time scaling for an 
A"-photon state will be rst(5', A^) oc . Hence, as pointed out by Steuernagel [30] . 
the exposure time will increase exponentially, and for realistic values of the number of 
detectors (pixels) and the efficiency of A^-photon absorbers, this will essentially rule out 
the practical use of quantum lithography. E.g., for a 5 x 5 pixel image the exposure 
time will multiply with 25 times going from a 2-photon state to a 3-photon state. For 
a more realistic image resolution of 100 x 100 pixels, the exposure time will increase 
by 10^ going from an A^-photon state to an A^ -|- 1 photon state. Therefore finding 
out which description is applicable is an important task to determine to what extent 
quantum lithography will be efficient. 

3. A conclusive experiment 

3.1. Description of the experiment 

For the experiment in figure [T]both theoretical explanations predict the same interference 
pattern when looking at two-photon absorption. Also for most other experiments 
performed so far (see for example [9] for double slits, [32l [33] for gratings instead of 
a double slit) either description could be used to explain the obtained results since 
only same position correlations were measured, and not the detection efficiency. To 
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Figure 2. Schematic experiment for measuring spatial correlations. See the text for 
more details. 



differentiate the theoretical models an experiment such as the one schematically depicted 
in figure |2] can be used. It is similar to the experiment in figure [1] in the sense that we 
ascertain that the state \'^2,si) is generated directly after the double slit. This can 
be achieved in the same manner as in the experiment previously described. However, 
instead of looking at same-position, two-photon correlation, we split the beam by an 
ordinary 50/50 beam splitter (BS) and detect the photons by two detectors that can 
have arbitrary positions relative the respective optical axis. We only look at the case 
where both detectors give a click, i.e., we are only looking at coincidences. We denote 
s as the position of one of the detectors and t as the position of the other detector. 
Subsequently the coincidence counts as a function of s and t are recorded. For such 
an experiment the two theoretical explanations predict different results. The detection 
probability of the two photons corresponding to ([5]) is 

of kdbs\ . . 

PB,sc{s,t) oc cos I — ^ I dst, (8) 

where 6st denotes the Kronecker delta. For ([6]), the detection probability becomes 

^st,sc(s,t) oc cos^ 

In figure [3] the detection coincidence rates are plotted as a function of s (horizontal) 
and t (vertical). (The lighter the colour, the higher is the coincidence rate.) Figure 
|3][a) illustrates the prediction of the model of Boto et al (j5]) and ([H]), where the photons 
retain spatial correlation. Figure W}^) illustrates the case of Steuernagel's description, 
corresponding to (E]) and ([9]), where the photons can propagate independently after the 
slit and hence allow an interference pattern even for the case when the two detectors are 
at different positions. The case where s = t, that is, along the figures' lower left to upper 
right diagonals, illustrate the results of most the experiments hitherto reported. For this 
case one can se that both theoretical models predict identical correlation patterns. 

3.2. Performing the experiment 

An experiment similar to that described in figure [2] was recently performed by Peeters et 
al ^28j in another context. The result of the measurement can be seen in figure 5(a) of 
their paper [28] and the figure is more or less identical to figure [3](b), indicating that the 
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Figure 3. Difference between tlie two theoretical explanations. In (a) the model 
of Boto et al is displayed. In (b) the model of Steuernagel is illustrated. The axes 
correspond to the positions s and t of the two detectors in figure [2] respectively. The 
distance between two peaks along the line s = i is half as long as it would be for the 
state in ([T]). The lighter the colour, the higher the coincidence rate is. The colour 
schemes in the two figures are only relative. One can see that for s = t the two 
theoretical models give the same pattern. 



two photons through the double sht are diffracted independently, but the two-photon 
character of the state will still be manifested in the coincidence correlation. 

Our understanding of the result is that since photons are non-interacting exchange 
particles, they will diffract and propagate independently in a linear medium. This 
means that in absence of some additional "guiding" , the photons prepared in the state 
|^2,si) will spread out according to normal diffraction theory, and it will be unlikely 
to detect the two photons at the same position far from the slit. This is also the 
case for more complicated states such as the inverse binomial states discussed in [11]. 
The entanglement will dictate the ensuing photon interference pattern, but each 
photon will in principle have its probability amplitude spread over the whole illuminated 
area, resulting in very low A^-photon, single detector, correlation. In Appendix B we 
show that this is also true for the slightly different setup discussed in [8]. In other 
words, Steuernagel seems to be correct in his analysis that both spatial and temporal 
concentration cannot be achieved at the same time at the required level for quantum 
lithography [34j- 

In principle, one could try to keep the spatial correlation by making the propagating 
media non-linear, and, e.g., form spatial solitons [35]. However, to manage this on a few 
photon level seems farfetched and not very realistic. 

Finally, one may ask if our conclusion is valid also for quantum lithography involving 
more than two modes [TOl [TTl [121 121]? The answer is yes. Invoking more modes will 
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allow the synthesis of more complicated patterns than a single slit interference pattern 
superimposed by a sinusoidal modulation, which is what can be achieved by a double slit. 
By adding more modes one can achieve additional superpositions of modes with different 
transverse wave vectors. Thus, one can make more involved lithographic patterns. 
However, if two photons emerging from one of two slits diffract independently, which is 
what the results of Peeters et al results indicate, then two (or more) photons will also 
diffract independently from any one of any number of slits or "holes" in a lithographic 
mask. Therefore, our conclusions about the exposure efficiency are applicable to any 
lithographic mask when combined with propagation through a linear media. 



4. Can spatial correlation be retained? 

One might wonder how our paper relates to other papers like |36j or J37j, where authors 
actually predict that A^-photon absorption proportional to J, where / is the normalised 
classical intensity of single-photon absorption, is theoretically possible? These two 
papers are also cited by Boto et al to motivate their argument that quantum lithography 
works with I instead of an dependence. For the results of these paper to be valid the 
authors postulate, without giving any specific mechanism, that the photons somehow 
are transported to exactly the same place, for example [36] states "We shall assume 
that, by focusing or otherwise, the photons of a correlated pair are brought at the same 
time onto a three-level atom..." . This, however, is not the situation in a typical quantum 
lithography scheme, since the photons are diffracted independently at a double slit (or 
more generally, in some lithographic mask) as our paper shows. In the proposal of 
Boto et al [8] and the experimental realisation of D'Angelo et al [9], which does not 
use a double slit, the photons are spreading independently as an analysis can show (see 



Appendix B ). Therefore, as long as light is propagating through linear media (such as 
air or vacuum), the results of [36] and [37] are not applicable in the existing quantum 
lithography schemes. 



5. Conclusion 



We showed that one can experimentally determine the scaling behaviour of the needed 
exposure time in quantum lithography as a function of the number of photons and 
the number of detectors (i.e., pixels) S. Such an experiment was recently reported, and 
it supports the view that photons, spatially entangled at a double slit to subsequently 
form a sub-wavelength pattern, will diffract independently. The photon state will still 
form the desired pattern (given that only two-, or in general, iV-photon coincidences 
are recorded). However, the experiment unfortunately also indicates that increasing the 
resolution through the number of entangled photons and by increasing the number of 
pixels one has to wait exponentially longer time to get the desired interference pattern. 
The explicit exposure time scaling is rst(>S', A^) oc 5*^. Therefore, quantum lithography 
involving a moderate number of pixels unfortunately seems to be impractical even for a 
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rather small number of photons. 
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Appendix A. Thorough analysis of the theoretical explanations 

Only assuming cw-pumping one can write the down-converted two-photon state as 

where the tilde sign~denotes Fourier transform and fc'/' (fcg") is the transverse component 
of the wave vector of the signal (idler) photon. E is the momentum representation of 
the pump profile in the crystal center plane and S is the phase matching profile. Their 
exact shape can be specified later. Transfering the above function from momentum to 
position space gives 

Vl^crystal«,4")0C / / d^i^ dfc^''Ep(fcr , fc'2'0S(fc;^ fe^'Oe-'^'^-^e'^^^''-^ , (A.2) 



where we have suppressed the Dirac notation with kets, as the information about 
subsequent interference lies in the probability amphtudes. The vector x'^' = {x"',y"'), 
i G {1,2} are the transverse coordinates in the crystal center plane and all the integrals 
here, and in the rest of the paper, are to be taken over the whole real axis. 

In the paraxial approximation, i.e., when the angle between signal and idler is not 
too large, one can write 

^,1,,{X'\,X",)0C J I d<d<Vl>„ystal«,0/^'(^"l,0/i'(a^"2,0,(A.3) 

where h' is the field propagator [38], which can be specified explicitly later on as well. 
The propagator h' allows any paraxial system including lenses and other components, 
but is especially simple for free-space propagation. The transverse coordinates just at 
the slit plane are denoted x"i. After the double slit we get 



^siit(fc"i,fc"2)oc / / d^^d^"2vI/,,,(^",,a."2)/«)/(^'2')e^'"^■""^e'^"-^"^(A.4) 




where k"e is the transverse wave vector of the signal {i = 1) or idler {£ = 2), just after 
the slit, but in the slit plane, and where we have assumed that the thickness of the slit 
is negligible. With a being the width of a slit, X < a < d, and d being the distance 
between the two slits, one can define 

_ J 1 for {d/2) - a/2 < \x';\ < {d/2) + a/2 
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i.e., the function is unity at the sht openings and otherwise. (The sht function / of 
course also depend on the width a, but we shall suppress that in the following as we 
shall be interested in relative, rather than absolute values of fields and intensities.) We 
have oriented the slits so that they are oriented along the y^'-direction and have their 
slit modulation structure in the x^'-direction. 

Finally we transfer the system to position space again and let it propagate to a 
detection plane so that we get 



*detect(a;i,a;2) a j J J y da;'ida;'2dfc"idfe"2*siit(A;"i, ^"2) 

X h{x^, x\)h{x^, x'2)e''^"'"'''Q~'^"^"''\ (A.6) 

where x'1^2 are the transverse coordinates directly after the slit and x-^ and X2 are 
the transverse coordinates of the detector positions for detecting signal and idler, 
respectively. Combining all the above formulae, one ends up with 

*detect(a3i,iC2) OC j dx i j dx 2 j dfc"i j dk"2 j dx"i j dx"2 

X j dx'i' I dx'^J dk'i' I dfe-4(fer,fc-)s(fcr,fc-) 

X h{x„ x\)h{x„ x'2)h\x'[, x";)h'{xl x'!;)f{x'[)f{xl) 

^ g-ifc'^■a,^g-ife''2•cc'2gifc'^■a;'^gifc''2•cc''2g-ifcf •<p-i^'^<_ 

Without loss of generality we can assume that the detectors are moved only 
perpendicular to the slits. It is then sufficient to evaluate all the integrals above only in 
the first X component of the vectors. Thus it suffices to write scalars instead of vectors, 
so in the paraxial approximation the field propagator can be written as 



h'{xlx-)^^^-^/^M-<) (A.8) 

and correspondingly for the other variables, where = z" — z'" is the distance of 
propagation. If /S.z is small then h'{x",x"') can be approximated by a delta-function 

5(x'i-x'n. 

Now we will make some assumptions. If the slits are very close to the crystal we 
have Az — >■ and therefore h'{x'(,x'(') — )■ 5{x'( — x'('). Assuming that the slits are 
narrower than any spatial structure of the two-photon field we can approximate 

« - d/2)S{x'^ - d/2) + S{x'; - d/2)5{x'i + d/2) 

+ 5{x'[ + d/2)5{x'!^ - d/2) + 5{x'[ + d/2)5{x'!^ + d/2), (A.9) 

where we have assumed that the two slits are identical. 

To proceed, the integration over k" can be simpfified as 

dfc^e^^^^'e-^^^^ ^ 6{xl - x'^). (A.IO) 



Combining all equations and carrying out the integrals leads to 

^(xi,X2) OC j (^K' J dk2Ep{k'l',k2 )E{k'l',k2) 
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X 



h{x^,d/2)h{x2, c//2)e-''=i''^/2g-ifc^"d/2 

+ -d/2)h{x^, rf/2)e''="''^/2g-ifc^"d/2 
+ h{x^, -d/2)h{x,, -d/2)e'Wei'=2"'i/2" 
Ciih{xi, d/2)h{x2, d/2) + ci-ih^x^, d/2)h{x2, —d/2) 
+ C2i/i(xi, -d/2)h{x2, d/2) + 022/^(^2, -d/2)h{x2, -d/2) (A. 11) 



where, e.g. 



Cll 



d*ryd.r4(*r.4")H(*r..;')e-'."-e-.^r.. 

dn / d'yEp(M,'j;)S(rf/2-u,d/2-t;). (A.12) 



In the following, since we discuss interference that require indistinguishability, 
we shall only consider the preparation of the symmetric state subspace at the slits. 
This requires symmetry of the functions E and S, implying that Cn = C22 and 
C12 = C21. Removing any global phase factor, and ignoring information about the 
absolute amplitudes, i.e. renormalising Cjj to Cij, such that |ci2p + jcup = 1, we can 
hence parameterise the amplitudes cn = C22 = e (2^9) sin (a/2), C12 = £21 = cos(q;/2). 
Apart from determining the probability amplitudes, the form of the wave function at 
the screen is independent of how the pump beam looks like and the particulars of the 
crystal, if one only imposes symmetry with the respect of the two slits. Of course, the 
particulars of the pump beam and the crystal determine the absolute magnitudes of Cn 
and C12. 

We want to stress that the requirement for symmetry (cn = C22 and C12 = C21) 
is automatically fulfilled if E and S are invariant under the simultaneously exchange 
of sign in both k'" and fcj". The pump beam profile E is usually a Gaussian beam 
depending only on the squared sum of the arguments, i.e. {k'" + ki^')'^, which makes it 
invariant under sign exchange, whereas the phase matching profile S under usual crystal 
orientation only depends on (fc'/' — fcg')^) making it invariant under sign exchange as well. 
However, it turns out that even under more advanced crystal orientations (the general 
case is given in equation (2) of [39]) the function S stays invariant under sign exchange 
of their arguments, although the dependence of the argument becomes more complex. 
Thus our symmetry requirement is almost naturally fulfilled for any crystal systems as 
long as the slits are symmetrical to the optical axis. 

Above, we assumed that the crystal is sufficiently close to the slit so that Az — t- 0. 
This assumption, however, can be dropped since any Az would transfer ^^crystaK^i', ^2") 
to some '^siit{k'(, ^2)5 and from there on, the analysis will be identical, ending up with 
the expression (lA.lip in this case too. This expression is in other words very general. 



The assumptions we cannot drop are that the slits are symmetrically excited by the 
two-photon field and that the field does not vary appreciably over the width of the slits. 



On the efficiency of quantum lithography 



13 



Inserting (lA.Sp into (lA.lip and considering a NOON state (a = vr) where both 
photons either take one sht or the other, substituting 2tt/X = k, and assuming a 
propagation distance R between the shts and the detector, we get 



The detection probabihty at the detectors hence becomes 

kd 



^NOON(a;i,a;2) oc cos 



(A.13) 



(A. 14) 



or, if Xi = X2, to PNOON(a;i) oc cos^ (^^i)- Equation (IA.14P is of the same form as 



([9]) (with bs = Xi and bt = X2), but now we arrived at it through a more thorough 
mathematical analysis. 

If we instead consider a state where the photons "chose" slit independently, so that 
Cii = C12 = C21 = C22, or a = 7r/2 and = 0, one obtains the coincidence probability 

P(Xi,X2) oc COS^ (l^^l) (y^^^^ ' ^^'^^^ 

or, if Xi = X2, to P{xi) oc cos"^ i^^i)- Comparing these expressions shows that the 
NOON-state leads to a doubling of the period, whereas allowing the photons to "chose" 
slit independently does not. 

Inserting ( ]A.8|) into ( lA.lip and varying the probability amplitude parameters a 



and (f, one obtains figure lAll These pictures correlate extremely well with figure 5 in 
the paper of Peeters et al |2Hj and show that the theory above provides an excellent 
description of their experiment. (Or rather, that Peeters et al are to be applauded for 
their ability to generate the various two-photon, two-slit states they aim for.) 

Note that these figures are not "ordinary" interference pictures, as they represent 
the correlations between two "point-detectors" along the same line, and not the intensity 
as a function of two-dimensional space. The idea of quantum lithography is illustrated in 
figure lAlT a). that describes NOON-state interference, flA.141) . and should be compared 
to figure [3l Single photon detection as in ([3]) can be seen as if one detector, let's say 
X2 is kept at a fixed position and only the other detector is moved. In that case one 
finds four and a half interference maxima along the xi axis. The two-photons absorption 
pattern is given by requiring that both detectors have the same coordinate, i.e., xi = X2- 
This corresponds to the diagonals in the figures, and going from the lower left to the 
upper right corner along the diagonal in (a) one counts nine maxima, i.e., a doubling of 
the number of fringes as compared to the single photon case. As a contrast, one does 
not see this period doubling for xi = X2 in (e), which corresponds to (1A.16I) . a state 
without spatial entanglement at the slits. This state will be produced, e.g., if the slits 
are located far from the crystal. In this case, the state after the slits becomes 



t + 6t)2 



vl>(xi, X2) = (|2, 0) + |1, 1) + |0, 2))/2 = ^ |0, 0) , (A.16) 

where (b^) is the creation operator of a photon through slit 1 (2). Hence, the 
transmission probabilities are those expected for two independent particles. As can 
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(b) 



(c) 






(d) 



(e) 



(f) 



Figure Al. Different plots of (|A.11| where each axis represents one detector 
coordinate, i.e., Xi and X2- A hghter colour indicate a higher intensity. The parameters 
are (a) a = tt, = 0, (b) a = 37r/2, ip = 0, (c) a = n/2, Lp = 7r/2, (d) a = 0, ip = 0, (e) 
a = tt/2, (f — 0, and (f) a — 7r/4, ip ~ 0. These calculated plots should be compared 
to the experimental results of Peeters et al , presented in figure 5 of their paper [2 8) . 



be seen from ( IA.15|) . the resulting two-photon interference pattern is the product of two 
single photon interference fringe patterns. 

The peculiar pattern in figure lAl( c) is the result of a superposition (|2,0) — 
i-\/2|l,l) + |0,2))/2 between a NOON state and the state in such a way that 

it cannot be written as the product of two independent photons. Figure lAl( d) 
describes the interference of the state and as expected, it does not show any 



same position (x^ 



X 



2>i 



two-photon interference but does show different position two- 



photon correlation fringes. However, what we ask the reader to take home from figure 
Al is that Steuernagel's description of quantum lithography describes the experiment of 
Peeters et al (correlations are predicted also where s ^ t), whereas the theory of Boto 
et al does not. 



Appendix B. Closer inspection of the original proposal of Boto et al 

We now assume a situation as in the proposal of Boto et al [8j and in the experiment 
of D'Angelo et al [9], i.e., two Gaussian beams meeting at a point z = 0, where the 
angles between the direction of propagation of the Gaussian beams and the z-axis is 
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±a and one of the beams is "above" the 2;-axis, whereas the other one is "below" so 
that the situation is symmetrical with respect to the 2;-axis. The two positions where 
the waists of the Gaussian beams are minimal are at a distance L away from the point 
X = y = z = 0, where we put the x-axis such that the two positions with minimum 
waists are at x = ±Lsin(a), y = 0, and z = — Lcos(a). The envelope function of a 
Gaussian beam at this minimum waist position can be written as 

-ikz I ■ ' 



^G(r-)=^oe-^'=^exp(-iA;-), (B.l^ 



where ipQ is a normalisation constant and 
1 1 . A 



q p nw"^ ' 



(B.2) 



where p the beam's phase radius of curvature, w the beam waist radius and r is the 
radial distance from the beam's optical axis. Moreover, at the mimimum waist position 
we have p = oo. 

The state for photons in a Gaussian beam is then given by 

\-^o{ru r-2, . . . , vn)) = ^pG{ri)tpG{r2) ■ ■ ■ MrN)al^al^ . . . 4^ |0, 0, . . . , 0) .(B.3) 

Here too, we shall suppress the Dirac notation in the following. The probability 
amplitude of the state after having propagated a distance L along its axis of propagation 
is given by 

^Lir'u r'„ rV) = e^^^^^ (^JY h{r\, r,)h{r',, r,) . . . h{r'^, r^) 

X 'ipG{ri)^pG{r2) ■ ■ ■ iJG{rN)dr^dr^ . . . drjs,, (B.4) 

where h{r' , r) is the field propagator. For Gaussian beams propagating a distance L 
these integrals can be solved and one can write 

v[/^(r'i, r'2, . . . , rV) = e-'''''^i^'a{r\Wa{r' 2) • • • V^^rV), (B.5) 

where [38] 

r 




V^gM = ^0 J^^exp (^-iA:^^-^J . (B.6) 

Coming back to the case of the proposal and the experiment, the detector is placed 
aX y = z = Q with x as the free variable and not with r' as in our description above. We 
therefore have to change the coordinates in ( IB. 61) . To do that, however, only photons 
arriving at a; = will have traveled the distance L, whereas photons arriving at any 
other position will have travelled the distance L =p AL(x), where the minus-sign is for 
the case of the Gaussian beam coming from above the z-axis and the plus-sign for the 
case, where the Gaussian beam comes from below the z-axis. One has the following 
geometrical relations: 

sin(a;) = AL/r', cos(a) = x/r\ tan(a;) = AL/x, (B.7) 
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where r' is the coordinate of r' which is in the plane spanned by x and z (the other 
coordinate of r' is parallel to y and will not change by our coordinate transformation). 
From now on we will for simplicity, but without loss of generality, set y = 0, since the 
detector is placed at y = and the coordinate systems of both two Gaussian beams 
have their y-axis oriented in the same direction. Equation (IB. 61) then transforms to 

[x cos(a)]^ 



exp 



-ik- 



(B.8) 



1 + [L =1= xsin(a)]/g \ 2[g + L =i= x sin(a)] ^ 

Since q in general is a complex number both the square root and the exponential 
function in (IB.SP are complex numbers. To be able to calculate any interference effects 
between the two Gaussian beams we will split the two functions explicitly into their real 
and imaginary parts. Introducing the dimensionless parameter /3 = X^L? / iji'^w'^) and 
looking only at the case when the minimum waist of the Gaussian beam is not too close 
to the detectors, i.e., L ^ xsin(a), we can simplify the above expressions to 



kx^ cos^ (a) ■ kx^ cos^ (a) sin (a) 
2L{l+/3-l)gT' 2L2(1+;9-1) .(B.9) 



The second and the last exponential-functions on the right hand side are the important 
terms here. The second term will be the Gaussian envelope and the last term will depend 
on whether a photon originated from the upper or the lower path. Having changed the 
coordinates of ipQ we have to be careful to also change the expression ^i-^^^^^) in (]B.4p 
and ( IB.SIl . since L there differs as well, depending on where the detector is placed. This 
leads to the replacement 



(B.IO) 

where the =F-sign is defined in the same manner as above. 

All the calculations so far have considered a single Gaussian beam. In the proposal 
of Boto et al , and in the experiment by D'Angelo et al , we have a NOON-state, i.e., 
a superposition of photons in an upper Gaussian beam and zero in a lower and vice 
versa. For such a case, the probability amphtude \l/ at the place of the detectors {z = 0) 
considering the results above can, after some algebra, be written as 



NOONla^l, X2, 



.Xn 



V2 



N 



X COS 



1 + py" 
I 

k sin(a) x 



g-iAfA;LgiAfarctan(v^)g- 



9 iV „ 

i — 1 

2L(l+^(-l) , 

N 



fccos*^(a) ^ 
2L(l+;3-l) 



v 



i=l 



k cos^(a) sin(a) ^ 

i=l 

2L2(l + /3-i) 



(B.ii: 



If we only look at the case of detecting all the A^ photons at the same place (i.e., 
Xi = X2 = ■ ■ ■ = Xn = x) we get the detection probability 



Pr 



NOONla;,x, . . .X) 



L(l + /3)i/4 



2N 



X COS 



e i(i+/3-'Ty 



k sin(a)A^x 



k cos (a) sin(a)A^x 
2L2(l + /3-i^ 



(B.12) 
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To see any interference effects we have to demand that the cosine-function oscillates 
faster than the exponential envelope function vanishes. This is fulfilled either for 

ksm(a)Nx > — — rr — B.13 

or for 

fccos^(a) sin(a;)A^a;^ k cos'^{a)y/PNx'^ m 

2L2(l + /3-i) ^ L{l + /3-^) ■ ^ ■ ^ 

The above relations are fulfilled if we have x < ^^^—r^^^ or x > 2nwy[X sin(a)], 
where we have used the definition of /3. For a > n/A one of the two conditions is always 
fulfilled. To get a feel for these terms, we can rewrite the cosine function in (IB. 121) as 

^x 27r cos^(a) sin(Q;)A^ /x\ 3 



cos^ 



27r sin(a)A^- 



(B.15) 



A 2(l + /3-i)LVA2 VA. 
and look at a typical example where we set a = 30°, A = 1 fim, L = 10 cm, w = 1 mm 
and X = 0.1 mm. For these parameters the pre-factor of a a;/A-term in the cosine 
function argument is 14 orders of magnitude higher than the pre-factor of the (x/A)^ 
term, so that the latter can safely be neglected in this case. 

Irrespective of which condition is fulfilled, the distance between two interference 
fringes scales with 1/A^ when going from 1 to photons and the width of the envelope 
function scales with 1 / a/ZV. This is the behaviour predicted by Boto et al . However, if 
we look again at ( IB. Ill) for the case N = 2 and with the assumption that L is sufficiently 
large, so that the x^-term in the cosine function can be neglected, we get 

p[xi,X2) = 2^^^ i(i+/3-M cos^ [A;sin(a)(xi + X2)] . (B.16) 

This function produces the same kind of pattern as in Steuernagel's model and shows 
that our somewhat simple description, not only describes the case of the double-slit 
experiment, but describes also the case of a NOON-state shared between two Gaussian 
beams, such as the setup-up of D'Angelo et al . In the derivation of (1B.16P we have 
assumed nothing but Gaussian beams and that the original state is in a NOON-state 
distributed over the two arms. 

We will now redo the previous analysis, but now assuming that all the photons are 
always produced at the same spatial coordinate. Instead of ( IB. 31) we use 

^0 (r-i, r2, . . . , vn) = ^^je-^'^^^e (-i^^) <^(^i - ^2)5(ri - rg) . . . 5(ri - r^) (B.17) 

as our input state probability amplitude, where t/'q is the new normalisation factor. This 
state is at the origin of the theory of Boto et al . The implication is that the first photon 
can essentially be found anywhere within the beam's Gaussian envelope, but the rest of 
the photons will subsequently be found at the same position as the first. 

In the same manner as before we can calculate the state after a distance L with 
help of the propagators. This leads to 



-N-l 



vl/^(r'i,r'2,...,rV) = ^^e-^^'^ ' ^ 



N + L/q\ XL 
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X exp 



— 1 



2L 



\ 



N 
i=l 



N 



N + L/q 



(B.18) 



which, by looking at the case where all the photons are detected at the same place, i.e.. 



r'2 



"'^ := (x, 0) and assuming the NOON-state leads to 

-N-l 



^NooN(a;) 



iNkL 



1 



iV2 + /3 



1/4 




XL 



2L(;3+JV2) 



2L(/3+JV^) 



xe 2L(/3+iv-) COS [A; (sin a) A^x] , (B.19) 

where we did the same replacements as in (]B.7p and (IB.lOp . setting ?/ = as before, 
and ignoring terms with higher orders than as before. 
This gives the A^-photon detection probability 



-PnoonI^;) 



1 

XL 



N-l 



cos^ [k (sina) A^x] .(B.20) 



This expression should be compared with (]B.12|) . where the x^-term is neglected. One 
sees immediately that both equations lead to the same distance between nearby peaks, 
given by the cosine-term. (This distance halving is what has been mainly been looked 
for in experiments so far.) Writing the equivalent to flB.16P when starting with fIB.lTp 
and flRTSj) leads to 



P(xi,X2) 



k cos^ (a)^/^(x^-\-X2)^ 
i{/3+4) 



cos'' [k sin(a) (xi + X2)] . (B.21) 



4 + /3AL 

Thus, in neither (IB.16P nor flB.2ip the photons are constrained to arrive at the same 
place as assumed by Boto et al . 

A difference, between (]B.12p and ( ]B.20|) is the exponential function. If one chooses 
L and w carefully so that one has /3 ^ 1, the argument of the exponential- function is 
proportional to oc w~'^ in (IB. 121) . but proportional to oc in (IB. 201) . Since 

the argument of the exponential-function determines the number of fringes one can see, 
one could, by changing w, determine which of the states one has produced, since in one 
case one would see more fringes, but in the other one less. This may be of interest for 
further experimental work in this field. 
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